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Abstract 
O ^ 

2 . We extensively develop a method of implementing mean-field calculations for de- 

formed nuclei, using the Gaussian expansion method (GEM). This GEM algorithm 
has the following advantages: (i) it can efficiently describe the energy-dependent 
^ . asymptotics of the wave functions at large r, (ii) it is applicable to various effective 

0> \ interactions including those with finite ranges, and (iii) the basis parameters are 

insensitive to nuclide, thereby many nuclei in wide mass range can be handled by a 
single set of bases. Superposing the spherical GEM bases with feasible truncation 
,j. , for the orbital angular momentum of the single-particle bases, we obtain deformed 

f^ ' single-particle wave-functions to reasonable precision. We apply the new algorithm 

OP ■ to the Hartree-Fock and the Hartree-Fock-Bogolyubov calculations of Mg nuclei 

with the Gogny interaction, by which neck structure of a deformed neutron halo is 
suggested for ^''Mg. 
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1 Introduction 

As experimental facilities supply fruitful data on nuclei far off the /^-stability, it has 
been recognized that theoretical approaches to nuclear structure should be renewed in 
some respects. One of the key ingredients is wave-function asymptotics at large r, which 
sometimes produces neutron (proton) halos in vicinity of the neutron (proton) drip line. 
It is also important to reinvestigate effective interactions, in connection to the magic 
numbers that are different from those near the /3-stability line. We developed a method 
for the spherical mean-field calculations [H [2] in which the Gaussian expansion method 
(GEM) [3] was applied. This method seems suitable to studying the shell structure from 
stable to drip-line nuclei, owing to its several advantages. 

A large number of nuclei have quadrupole deformation. The deformation plays a sig- 
nificant role in unstable nuclei as well. For instance, the deformation can be relevant to 
the disappearance of the N = 20 magic number on the so-called 'island of inversion' [1], 
whereas there remain arguments for individual nuclei; e.g. spherical description has also 
been proposed for '^^Mg [5j. It is desired to implement calculations with deformation, 
reproducing the wave-function asymptotics and handling a wide variety of effective inter- 
actions simultaneously. 

In this paper we propose a new method for deformed mean-field calculations by ap- 
plying the GEM. Taking advantage of the flexibility in describing the radial degrees of 
freedom, we adopt a set of the spherical GEM bases to represent deformed single-particle 
(s.p.) wave functions, with truncation for the orbital angular momentum i. The new 
method is tested in the Hartree-Fock (HF) and the Hartree-Fock-Bogolyubov (HFB) cal- 
culations for Mg nuclei with the Gogny interaction, and the results are compared with 
those in literatures. From the present calculation, an interesting feature is suggested for 
the drip-line nucleus ^'^Mg. 

2 Single-particle bases 

In this paper we assume the nuclear mean fields to be axially symmetric and to conserve 
the parity. The z axis is taken to be the symmetry axis. The method can immediately be 
extended to general cases with no symmetry assumptions on the one-body fields, apart 
from an additional constraint on the center-of-mass (cm.) position. 

2.1 GEM bases 

We represent the s.p. wave functions by superposing the spherical Gaussian bases, which 
have the following form: 

Ruij{r) = A/;^j/exp(-z/r2). ^-^^j 

Here Y^^'{r) expresses the spherical harmonics and Xo- the spin wave function. We drop the 
isospin index without confusion. The range parameter of the Gaussian basis z/ is a complex 
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number in general [6]; i^ = z/r + iz/j {uj- > 0). Via the imaginary part oscillating behavior of 
the s.p. wave functions can be expressed efficiently [2]. Formulae for calculating the one- 
and two-body matrix elements that are required in the HF and the HFB calculations, as 
well as the constant A/'^^j, are given in Refs. [H [2]. The s.p. wave functions under the 
axially deformed mean field are represented as 

1pnnm{r) = Y. cS]m <^v£jm(r) , (2) 

where the subscript vr on the Ihs stands for the parity. The sum of i and j on the rhs 
runs over all possible values satisfying it = {—Y, j = i i: 1/2 and j > \m\, in principle. 

In the GEM we usually take z/'s belonging to a geometric progression. In Ref. [2], 
we found that a certain combination of the real- and complex-range Gaussian bases is 
suitable for nuclear mean-field calculations. In all the following calculations, we take the 

basis-set of 

( u, = (n = 0,l,---,5) 

-. = -ofe-^", il^i- (n = 0,l,2) ' (3) 

with z/q = (2.40 fm)~^ and h = 1.25, irrespective of {i,j). Namely, 12 bases are employed 
for each {i,j)', 6 bases have real u and the other 6 have complex u. This set is quite 
similar to Set C in Ref. [2]. 

2.2 Adaptability to wave functions with various size 

An appropriate set of the GEM bases is capable of describing wave functions with various 
size. This feature is desirable for self-consistent mean-field calculations with deformation, 
because in deformed nuclei the density distribution depends on the direction, and degree 
of the deformation is not known in advance. We here show adaptability of the GEM 
with respect to size of nuclei, by presenting results of the spherical Hartree-Fock (HF) 
calculations. 

Size of spherical nuclei, which is typically represented by the rms radii, depends on 
the mass number A, apart from exotic structure such as neutron halos near the drip line. 
Many mean-field calculations have been implemented by using the harmonic oscillator 
(HO) bases, particularly when the effective interaction has finite ranges. In the mean- 
field calculations with the HO bases, the length parameter of the bases 6^ = l/vMcJ 
depends on A. For stable nuclei, 6^ is almost proportional to A^^^, as a; ~ 41.2 A~^^^ MeV, 
although bi^ is often adjusted for individual nuclides so as to minimize their energy E. 
In contrast, since the GEM basis-set contains Gaussians of various ranges, even a single 
set can describe many nuclei to good precision. Binding energies and rms matter radii 
calculated with the Gogny DIS interaction |7j are tabulated in Table [H for the doubly- 
magic nuclei ^"^O, ^^Q, ^^Ca, ^^Ca, ^°Zt and ^o^Pb. The values obtained from the GEM 
basis-set of Eq. ([3]) are compared with those from the A-dependent HO basis-set. The 
Coulomb interaction between protons is handled exactly [1], and the cm. motion is fully 
removed from the effective Hamiltonian before variation. The influence of the cm. motion 



Table 1: Binding energies —E (MeV) and rms matter radii J (r^) (fm) of ^^O, ^^O, ^°Ca, 
^^Ca, ^°Zr and ^osp]^ ^^ ^^le HF calculations with the DIS interaction. Results of the 
GEM bases and the HO bases are compared. 



Nuclide 




HO 


GEM 


WQ 


-E 


129.638 


129.520 




\/{r') 


2.603 


2.606 


240 


-E 


168.573 


168.598 




\/{r') 


2.998 


3.014 


40Ca 


-E 


344.470 


344.570 




\Kr') 


3.370 


3.370 


48Ca 


-E 


416.567 


416.764 




\/{r') 


3.514 


3.513 


902r 


-E 


785.126 


785.928 




\/{r') 


4.238 


4.235 


208pb 


-E 


1638.094 


1639.047 




^(r^) 


5.517 


5.513 



on the rms matter radii is treated in a similar manner [Hj • In the calculations using the HO 
bases, the 6^ parameter of the bases is determined from to = 41.2 A~^^^ MeV, and all the 
bases up to A'osc = 15 are included, where A^osc is number of the oscillator quanta. Because 
of the variational nature of the HF theory, the lower energy indicates the more reliable 
result for individual nuclei. In this regard, the A-independent set of the GEM bases gives 
no worse, even slightly better results than the A-dependent HO basis-set except for ^^O. 
To show further the adaptability of the GEM basis-set to nuclear size, the calculated 
density distributions of the six nuclei are illustrated in Fig. [1] 

In the results of Table [T|, we use larger number of bases in the GEM calculation than 
in the HO calculation. Although one might think that this is not fair comparison, it 
is not technically easy to increase the number of the HO bases, because of round-off 
errors. It is rather another advantage of the GEM that we can avoid harmful round-off 
errors even with a large number of bases. Furthermore, it should be stressed that the 
GEM basis-set used here is nucleus-independent, demonstrating adaptability to nuclear 
size. This yields an additional advantage in systematic calculations. In the self-consistent 
mean-field calculations such as HF and HFB with finite-range interactions, computation 
of the two-body interaction matrix elements is the most time-consuming part. Since the 
basis-set of Subsec. 12. H is applicable to many nuclei, we do not have to repeat computation 
of the interaction matrix elements once they are stored. 
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Figure 1: Density distribution of ^^O (black line), ^^O (brown line), ^°Ca (blue line), '^^Ca 
(green line), ^°Zr (purple line) and ^osp]^ ^^.^^ line), calculated by the GEM basis-set with 
the DIS interaction. 

2.3 Test on truncation with respect to i 

As viewed in Eq. (|2]), there is mixing of i and j in the s.p. wave functions under the 
deformed mean field. In particular, the ^-mixing is relevant to the deformation. The 
mixing of j (= i ± 1/2) is driven by the spin-orbit coupling. In deformed mean-field 
calculations with the spherical bases, we necessarily truncate the s.p. bases with respect 
to i (and j), by restricting the sum over i in Eq. ([2]) as £ < £max- It is an unavoidable 
question how large ^max is needed for practical calculations. 

To answer this question, we consider an axially symmetric HO potential, giving a s.p. 
Hamiltonian of 
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where cUq 



{2ujI + cu2)/3 and S, 



def 



{uj_ — uj1)/2uJq. The exact eigenvalues of this 
Hamiltonian are obviously e{nx_nz) = ui±{n_i + 1) -^(^^(n^ -|- i) (nj_, n^ = 0, 1, ■ ■ ■ ). On the 
other hand, we obtain approximate eigenvalues by diagonalizing this Hamiltonian with 
the spherical GEM bases, by setting a certain value of £max- Quality of the £max truncation 
can be assessed by comparing the approximate solutions to the exact ones. 

In Fig. [2], we plot errors of the energy eigenvalues obtained by the GEM, Ae{n±nz) = 
^GEuin^n;,) — e exact {n±nz), as a function of Sdd, taking £niax = 7 with the basis-set of 
Eq. i^. Although we set {ulu^y/^ = 41.2 A~^/^ MeV with A = 24 at each 6dei, the errors 
are insensitive to uo if we measure Ae in unit of uq. We have confirmed that the results 
presented in Fig. [2] are almost unchanged even for A = 208. The eigenvalues belonging to 
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Figure 2: Errors of the GEM solutions in energy eigenvalues of the anisotropic harmonic 
oscillator of Eq. (j4]), whose ranges are represented by the hatched areas for each A'osc 
value. Upper panel: Errors for the vr = + levels. Red, blue and green areas show bounds 
of errors for the A'osc = 0,2 and 4 levels, respectively. Lower panel: Errors for the vr = — 
levels. Red and blue areas represent bounds of errors for the A^osc = 1 and 3 levels. 



equal A^osc(= n± + n^) have errors of the same orders of magnitude (except at ^dcf = 0). 
This is because the levels belonging to equal A'osc significantly mix one another owing to 
their close energies, while admixture of the iVosc ± 2 levels occurs only perturbatively. We 
therefore classify the s.p. levels by iVosc and show their errors as a function of (5def by the 
hatched areas, which are edged by the minimum and maximum errors for each iVosc- 

At (5dcf = {i.e. ujx_ = LJz), the energy eigenvalues become e{n±n;^) = ujq{Nosc + §)• 
The errors at 6dci = exclusively come from the choice of the radial parameters in Eq. ([3]) , 
irrelevant to the ^-mixing. We have much smaller errors at Sdd = than those at finite 
Sdei- Because the radial parts are well described, the errors at (5dcf j^ predominantly 
depend on how well the ^-mixing is taken into account. Irregularities at 5dcf > 0.5 as 
well as at ^dcf = occur due to crossing of the levels that have equal vr and m quantum 
numbers. 

Whereas it depends on the situation what precision is required in the calculation, we 
here assume the criterion to be Ae < 0.01 cuq? which corresponds to Ae < 0.1 MeV for 
nuclei with uq ~ 10 MeV. With the £max = 7 truncation, all the A'osc < 3 levels satisfy 
this criterion at —0.5 < 6dci ^ +0.35. This region of 6dcf covers deformation of most 
nuclei, as long as their low-energy states are concerned. Since Ae exceeds 0.01 cuq even 
at Sdcf ~ 0.1, we cannot expect that the £max = 7 truncation works well for the iVosc > 4 
levels. It should be noticed that the crossing between the highest level of iVosc = 3 and 
the lowest level of iVosc = 5 takes place at 6dcf = +0.41 and —0.56. These crossing points 
are close to the edges of the region where the above criterion is fulfilled. Beyond these 6dei 
values some of the A^'osc = 5 levels, for which we need higher £max, are occupied before all 
the iVosc = 3 levels are, although the crossing levels do not necessarily mix since they may 
have different quantum numbers m. Similar crossing between the TVqsc = 4 and 6 levels 
occurs at Sdei ~ +0.34 and —0.44. Thus the crossing points between A^osc and A'osc + 2 
levels are not very sensitive to A'osc- Although the spin-orbit splitting influences the level 
crossing points, the current argument will hold in more realistic cases to a good degree, 
because the spin-orbit splitting is much smaller than 2ci;o. Thus it will be reasonable to 
state that we should take ^max > N^^^ + 4 in deformed mean-field calculations, where A^^^ 
is defined to be the highest oscillator quantum number for the s.p. levels under interest. If 
we take £max > A''^^ + 6, we can expect better precision by about two orders of magnitude; 
the error on the level corresponding to N^^^ will be Ae^ < lO'^c^o for low energy states. 

While we have discussed errors on the s.p. energies, mean-field calculations are carried 
out by minimizing the total energy of each nucleus. Because the errors grow for the larger 
\6dei\, calculations with the spherical bases tend to yield higher energy for deformed states 
than for spherical states. We put A^^^^ to be the value for the Fermi level, which can be 
estimated as A^^^^ ^ (374/2)^/"^ |9]. The number of nucleons occupying the A^^^, levels 
should be 2{N^^^Y ^ {3\^A)'^^^ . Dominated by the errors on the A^^^^ levels, error on 
the total energy due to the £max truncation is estimated to be AE ^ {3\/2AY^^ Ae^. 
Inserting Ae^ ~ O.OIuq and luq ^ 40 MeV, we obtain AE «i A^^^ MeV for well-deformed 
configurations, while this error is absent in spherical configurations. We should keep this 
point in mind when comparing energies of a spherical minimum and a deformed one. 



3 Implementation of self-consistent mean-field calcu- 
lations with deformation 

We have discussed that the spherical GEM bases shown in Subsec. 12.11 is promising in 
describing deformed nuclei. We next present how we implement self-consistent HF or 
HFB calculations using the GEM bases. 

One of the advantages of the GEM algorithm in Refs. [11 [2j was that finite-range 
effective interactions are tractable [8] by computing two-body interaction matrix elements 
of the interactions, even for the LS and the tensor channels. This advantage is maintained 
in the deformed cases. While in the spherical cases the number of necessary matrix 
elements is reduced to great extent due to the symmetry of the one-body fields, we need 
all the non-vanishing two-body matrix elements in the deformed mean-field calculations. 
Since we adopt the spherical GEM bases, the matrix elements can straightforwardly be 
computed according to the formulae shown in Refs. [HE]- It costs 4.3GB of memory 
or disk to store the matrix elements with i < 7 in double precision. We handle the 
Coulomb interaction in a similar manner, which needs additional 2.2 GB. Although it 
is time-consuming task to compute them, those matrix elements are useful for various 
calculations once we store them; not only for HF and HFB calculations of many nuclei, but 
also for calculations in the random-phase approximation, as will be discussed elsewhere. 

The density-dependent interaction should be renewed at each iteration. This is not a 
difficult task as long as this part of the interaction has a contact form as in Ref . [8] . 

Folding the stored matrix elements with given s.p. wave functions, we construct the 
s.p. Hamiltonian that preserves the (vr, m) quantum numbers. The s.p. wave functions are 
then obtained by solving the HF or HFB equation as a generalized eigenvalue problem, for 
each (vr, TTi). Starting from appropriate initial values, we repeat this procedure iteratively 
until convergence. Alternatively, we can apply the gradient method to obtain the energy 
minimum by using the s.p. Hamiltonian. 

The present algorithm keeps the advantages of the GEM in the spherical mean-field 
calculations. We here list them again: (i) it is efficient in describing the energy- dependent 
asymptotics of s.p. wave functions at large r, (ii) we can handle various effective inter- 
actions, including those having non-locality, and (iii) the basis parameters are insensitive 
to nuclide, thereby a single-set of bases is applicable to wide mass range of nuclei. 

4 Numerical examples for magnesium isotopes 

We now apply the present method of the deformed HF and HFB calculations to actual 
nuclei. We take even-A^ magnesium isotopes as examples, which contain a well-deformed 
stable nucleus ^^Mg, several nuclei on the island of inversion such as ^^Mg and "^^Mg, 
and '^^Mg for which the N = 28 magicity has been predicted to disappear in several 
calculations. We exactly treat the Coulomb energy [1], and both the one- and two-body 
terms of the cm. part are subtracted from the Hamiltonian before variation, unless 
mentioned explicitly. The GEM bases of Subsec. 12.11 are employed with the fmax = 7 



truncation. This fulfills £max > N^^^ + 4 for the ground states of all the Mg nuclei to be 
presented, since the last neutron occupies an N^sc = 3 orbit in the neutron-rich Mg nuclei. 

4.1 Comparison with previous calculations 

We examine how well the present method can describe the deformed Mg nuclei. We first 
show results of the HF calculation with the Dl parameter-set [10] of the Gogny interaction. 

In Ref. IH], an axial HF calculation using the HO bases has been implemented for 
many Mg isotopes. Although the basis-set consists only of A'osc < 4, the oscillator length 
6^ is optimized except for '^^Mg and ^^Mg. In Table[2]we compare the present GEM results 
for the Mg isotopes with those in Ref. [11]. For ^"^^^^^^^^Mg, binding energies computed 
in the antisymmetrized molecular dynamics (AMD) [12] are reported. The AMD is a 
powerful tool to study structure of light to medium-mass nuclei. In the AMD approach 
the total wave function is represented by a Slater determinant of Gaussian wave packets 
of constituent nucleons. While ranges of the Gaussians are taken to be equal for all 
nucleons in the nucleus, the nucleons have different central positions from one another 
and the positions are optimized without assuming the axial symmetry. Whereas the 
AMD has been extended by incorporating the projections and superposing many Slater 
determinants, we here use the results of its simplest version for comparison, because it is 
analogous to the HF approximation. 

With respect to the spurious cm. motion, it is popular in the mean-field calculations 
so far to remove only the one-body term before variation. One of the present GEM results 
(denoted by GEMi) are obtained by this prescription. While infiuence of the spurious 
cm. motion can be fully removed in the AMD, only the one-body term is subtracted 
before variation in the AMD results shown in Table El Therefore the GEMi values should 
be compared to the HO and the AMD ones. We additionally present the values in which 
the two-body term of the cm. Hamiltonian is also subtracted before variation in the 
GEM calculation, and denote them by GEM2. 

Since the HF approximation holds variational nature, the lower ground-state energy 
indicates the more reliable result {i.e. the closer to the true minimum). In this respect 
the present GEM bases give more favorable results than the HO and the AMD bases. The 
present method gives 7 — 8 MeV lower energies than the AMD in '^^''^'^''^^Mg. Compared to 
the HO calculation in Ref. [11], the energy gain in the present calculation grows as going 
to the neutron-rich region, ~ 2 MeV in ^^Mg to ~ 10 MeV in ^^Mg. It may be due to 
spatially broad distribution of density in the neutron-rich region, which is well described 
by the GEM but is hard to be reproduced with the HO bases. It should be noticed that 
38]VIg is unbound in the result of Ref. [11] because it has higher energy than '^^Mg, while it 
is bound in the present calculation. This exemplifies importance of numerical algorithm 
that appropriately handles spatial extension of wave functions when investigating the 
neutron drip line. 

The two-body term of the cm. Hamiltonian affects the energies to sizable amount, 
as clearly viewed in Table [2l Difference between the GEMi and GEM2 energies slightly 
grows for increasing A. 



Table 2: Comparison of the axial HF results for binding energies —E (MeV) and intrinsic 
mass quadrupole moments Qq (fm^) of the Mg isotopes. The present GEM results are 
compared with those of the HO bases [11] and the AMD calculation [12] . The one-body 
cm. energy is subtracted in GEMi as well as in the HO and the AMD results, while the 
two-body cm. energy is also removed in GEM2. The Dl interaction is adopted for all 
the calculations. 



Nuclide 




HO 


AMD 


GEMi 


GEM2 


22Mg 


-E 


— 


161.26 


169.03 


162.44 




Qo 


— 


— 


100.84 


104.53 


24Mg 


-E 


193.90 


188.33 


195.73 


188.74 




Qo 


113.95 


— 


114.39 


118.11 


26Mg 


-E 


211.11 


— 


212.99 


205.29 




Qo 


76.51 


— 


78.20 


81.09 


28Mg 


-E 


227.82 


223.07 


230.14 


222.45 




Qo 


101.37 


— 


101.41 


105.06 


30Mg 


-E 


237.69 


— 


240.34 


232.36 




Qo 


76.10 


— 


75.07 


77.29 


32Mg 


-E 


245.98 


— 


249.05 


240.97 




Qo 


34.27 


— 


37.13 


38.33 


34Mg 


-E 


249.94 


— 


254.77 


246.43 




Qo 


111.70 


— 


128.36 


131.92 


36Mg 


-E 


253.79 


— 


260.49 


251.84 




Qo 


146.17 


— 


176.50 


180.81 


38Mg 


-E 


251.69 


— 


261.69 


252.62 




Qo 


168.65 


— 


175.06 


179.47 



Table 3: Binding energies -E (MeV) of 30,32,34j^g ^^ ^^^ gpg calculations with the DIS 
interaction. The GEM results are compared with those of the HO bases in Ref. J13j . 

^Nuchde HO GEM Exp. ^ 



30Mg 239.30 239.47 241.63 
32Mg 248.22 248.30 249.69 



34 



Mg 252.82 254.02 256.59 



We next turn to the HFB calculations. We compare the GEM results for ^''"^^Mg 
with those obtained from the HO bases truncated by A^osc < 10 [I3], in Table [3l The 
DIS parameter-set [7j of the Gogny interaction is adopted, and the cm. Hamiltonian is 
fully subtracted before iteration in both calculations. In Ref. [13] , the spherical HO bases 
are used, so that the angular-momentum projection could be carried out afterward. The 
exchange term of the Coulomb energy is handled in the Slater approximation in Ref. [13] . 

We find that the present method gives lower energies in all of '^^"^^Mg than the HO 
calculation in Ref. [13]. If we take differences in computational procedure irrelevant to 
the basis-sets {e.g. the treatment of the Coulomb exchange term) into consideration, it is 
fair to say that the present results are no worse than the HO ones in Ref. [13] from the 
variational viewpoint. We obtain lower energy in the present calculation than in Ref. [I3] 
by greater than 1 MeV for the unstable nucleus ^'^Mg, while by less than 0.2 MeV for 
30,32]y[g^ This suggests that the spatial extention is more important than the mixing of 
the i > 7 components particularly in ^^Mg. Notice that the HO basis-set in Ref. [13] 
contains up to ^ = 10. 

4.2 Neck structure of neutron halo in ^^Mg 

As we have pointed out, the GEM can describe the wave- function asymptotics at large 
r to reasonable precision. Taking this advantage, we investigate density distribution of 
axially deformed drip-line nuclei, within the mean-field approximation. 

We here discuss in the HF framework. The s.p. wave function ipniTm{T^) in Eq. ([2]) can 
be decomposed by a sum of spherical wave functions, 

ll^nnm (r) = ^ cJaL vtm (l") • (5) 

Here <^^^m(^) — ^e^Li''") ^me i^) Xa with m^ + a/2 = m, which is defined by the (icr)- 
projection on ipnnm{T') besides a constant factor. The nuclear force becomes negligible 
at sufficiently large r. There the coordinate-represented HF Hamiltonian for neutrons 
becomes approximately spherically symmetric and the components having different i de- 
couple to one another [H]. The spin-orbit coupling, which is an effect of the nuclear 
force, also becomes negligible, and thereby the a (and mi) value is frozen. After each 
(ia) component propagates over a certain region of r, ^„7rm(r) will be dominated by the 
'^la-mi^) component of the lowest possible i {i.e. i = \m\ ± 1/2, with the sign fixed by 
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the parity) due to the centrifugal barrier [TU [15] , although the degree of the dominance 
depends on characters of the deformed orbit. The asymptotic form of R^^lnir) at large r 
is e~''"'''"'"/r with rjmrm = \/2M|e„^m| i^mrm is the s.p. energy). Hence the neutron density 
distribution asymptotically behaves as 



Pnirj 



^p|#(r)p^^p^5:ic(-) 



£p a rrip I 






(6) 



where the subscript F represents the highest occupied level, ip is determined by ip = 
Imp I ± 1/2 and (— )^^ = vrp, and gp stands for the occupation number on the level F. We 
hereafter assume mp > without loss of generality, owing to the reflection symmetry. 

In connection to halos, the £p = and 1 cases are particularly interesting. For ip = 0, 
which implicates Tip = + and mp = 1/2, Eq. (Q is reduced to 



p-'iripr 
f \ ^. |-("f) 1 2 

P„(r)~a^^:^|c^+;/2l 



Y^\r] 



,2 l'^0 + l/2l ' 



All 



r^ 



(7) 



giving isotropic density in the asymptotic region. Hence halos formed by vr = + last 
nucleon are spherically symmetric, as long as the i > 2 components can be neglected. 
For ip = 1 which presumes iip = — , there are two possibilities mp = 1/2 and 3/2. In the 
mp = 3/2 case, Eq. (^ becomes 
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This asymptotic component is obviously deformed, having vanishing contribution in the 
z direction. In the m,p = 1/2 case, we obtain 
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The asymptotic component of Eq. Qj can also be deformed, and its degree depends on 



'1 + 1/21 



and |cj^"^{/2l- Thus, within the HF approximation, a deformed halo is expected 
in axially deformed drip-line nuclei with the last nucleon occupying a {7T,m) = (— , 1/2) 
or (-,3/2) orbit. 

Connected to the Nilsson model, the asymptotic quantum numbers have widely been 
used in description of well-deformed nuclei, comprising m,i (denoted by A in the Nilsson 
model). Therefore a becomes an approximate quantum number. If this is the case, one of 
I ^1^+1/2 1 ^^ much greater than the other in Eq. (Q. Then the asymptotic behavior of Pn(r) 



11 



depends largely on direction in the {nF,mp) = (— , 1/2) case. In the extreme case where 
one of clXi/2 vanishes, the damping factor e~^'^^/r^ is missing in a certain direction; if 
'^1-1/2 ~ 0' Pn damps according to the e'"^"^^ /r'^ factor in the z direction, but damps 
faster in the x and y directions. The absence of the e'"^^^ /r"^ asymptotic component is 
significant only in the vicinity of the xy plane. Hence a neck structure of a halo is expected 
for drip-line nuclei with [np, mp) = (— , 1/2) that is dominated by the rrii = component. 
In contrast, the density in nuclei with {7ip,mp) = (—,3/2) or (—,1/2) dominated by 
m£ = 1 may have dips on the z axis, forming a halo with an apple-like shape. 

In the Mg isotopes, the neutron-drip line has been predicted to lie at A^ = 28 in several 
mean-field calculations so far [161 Ej- ^^ ^ recent experiment '*°Mg has been confirmed 
to be bound [18]. We obtain the drip line at A^ = 28, as reported in Ref. [IT]; "^"Mg 
is bound while ^^Mg has higher energy than '^'^Mg. In the present HF result ^"Mg has 



a prolate shape at its energy minimum, with E = — 260.7 MeV, J {r^) = 3.63 fm and 
Qq = 227 fm . It is noted that this energy is appreciably lower than in the A'^osc < 17 HO 
calculation [17]. This is probably connected to the spatial extension of the wave function 
shown below. The last two neutrons, which have ep = —1.14 MeV, occupy a (vr, m) = 
(— , 1/2) orbit that has large portion of the £ = 1 component by 68%. Moreover, we have 
l4"f i/2p/(l^iTi/2l^ + 1*^1-1/2^) ~ 0.95 for this orbit. Applying the above argument, it is 
likely that Pn(r) of ^°Mg has a halo oriented to the z direction, but this halo component 
is almost missing on the xy plane, leading to neck structure of the neutron halo. 

In Fig. [3l pr{x = y = 0,z) and Pr^x^y = z = Q) (r = p,n) of ^°Mg are depicted. 
Although Pp{z) is elongated in comparison to Pp{x) to certain extent because of the 
prolate deformation, they damp with nearly equal slope at r > 5fm. We find halo 
structure in Pn{z) at 2 > 8fm with the asymptotic behavior that is consistent with ep 
to good approximation. It is remarked that the component with the same asymptotics is 
highly suppressed in the x direction. 

To observe the shape of "^"Mg at low density, equi-density lines on the x2;-plane are 
drawn in Fig. |H For comparison, similar plots for ^^Mg and ^^Mg are given as well. All the 
three Mg nuclei have prolate deformation in the present HF calculation. The equi-density 
lines distribute with almost equal intervals from 10~^ to 10~^ fm~^, indicating exponential 
decrease of Pr- With the exception of pn in ^''Mg, the size of the intervals does not depend 
on direction. Although the p^ = 10~^fm~^ line in ^^Mg is slightly constricted in the x 
direction, it could be influenced by deviation from the exponential asymptotics due to 
numerical errors. The contour plot of p„ in ^°Mg is intriguing. The large intervals of the 
equi-density lines correspond to the slow damping of p„ due to the small ep. However, the 
intervals are comparable to those in ^^Mg in the x direction. We view the neck structure 
in the low p„ region. 

Because it mixes the m£ = and 1 components, the pair correlation tends to restore 
the spherical symmetry in the density distribution in the asymptotic region. Depending 
on degree of the mixing, the neck structure may survive or disappear. However, in the 
present HFB calculation with the DIS interaction, we do not find a superfiuid solution 
having lower energy than the HF minimum in ^°Mg. Thus the density distribution of 
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Figure 3: Proton (in red color) and neutron (in blue color) densities of ^''Mg in the x 
(dashed lines) and the z (solid lines) directions, calculated with the DIS interaction. 

■^^Mg shown in Figs. [3] and S] is not altered. 

If the a > 2 components remain sizable in the asymptotic region, deformed halos 
(including neck structure) with other {11,771) may be present. In Ref. [19], possibility of 
deformed halos depending on the direction was pointed out for ^^Be and ^^C, based on 
the Skyrme HFB calculations. Those halos are formed by a (vr, m) = (+, 1/2) orbit. 

There remain problems with respect to deformed halos: (a) how correlations beyond 
HFB (including restoration of the rotational symmetry) affect them, and (b) whether 
they are detectable in experiments. Although both questions are very important, it is not 
easy to give satisfactory answers at this moment and we leave these problems to future 
studies. We here emphasize that such exotic structure of nuclear halos can be investigated 
only via numerical methods that are capable of describing the wave-function asymptotics 
appropriately. 

5 Summary 

We extensively develop a new method of implementing the Hartree-Fock (HF) and the 
Hartree-Fock-Bogolyubov (HFB) calculations of nuclei with deformation, applying the 
Gaussian expansion method (GEM). Owing to the adaptability in describing radial de- 
grees of freedom, which is confirmed by the HF calculations from ^^O to ^'^^Pb, we adopt 
the spherical GEM bases. We argue how large i should be taken into account, by compar- 
ing the numerical solutions obtained from the spherical GEM bases with the analytic ones 
in the axially deformed harmonic oscillator. The present method maintains three notable 
advantages of the GEM algorithm for the mean-field calculations: (i) we can efficiently de- 
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Figure 4: Contour plot of p^. (r = p, n) on the xz-plane for ^^Mg, ^^Mg and ^°Mg, 
calculated with the DIS interaction. Positions of Pr = 0.05, 10~^, 10~^, lO""^, 10~^ and 
lO^^fm" are presented. 
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scribe the energy-dependent asymptotics of single-particle (s.p.) wave functions at large 
r, (ii) we can handle various effective interactions, including those having non-locality, 
and (iii) a single-set of bases is applicable to wide mass range of nuclei and therefore is 
suitable to systematic calculations. 

The present method is applied to magnesium nuclei with the Gogny force. We show 
that the present results are no worse than those in the literatures, from the variational 
viewpoint. Compared with the conventional method using the harmonic oscillator bases, 
the present method is suitable particularly to nuclei far from the /3 stability. For ^°Mg, 
we suggest neck structure of a neutron halo, which arises due to the asymptotic behavior 
depending on the direction. Such a possibility can be argued only via methods describing 
the asymptotic form of the s.p. wave functions appropriately. 

This work is financially supported in part as Grant-in-Aid for Scientific Research (C), 
No. 19540262, by Japan Society for the Promotion of Science. Numerical calculations are 
performed on HITAC SRllOOO at Institute of Media and Information Technology, Chiba 
University, and at Information Technology Center, University of Tokyo. 
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